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^ '. Abstract 



> 

\^ ' The roughening of interfaces moving in inhomogeneous media is investi- 

^O ■ gated by numerical integration of the phenomenological stochastic differential 

Q^ I equation proposed by Kardar, Parisi, and Zhang [Phys. Rev. Lett. 56, 889, 

-4— > ' 

d ' (1986)] with quenched noise (QKPZ). We express the evolution equations for 

the mean height and the roughness into two contributions: the local and the 

lateral one. We compare this two contributions with the ones obtained for 

two directed percolation deppining models (DPD): the Tang and Leschhorn 



X ■ model [Phys. Rev A 45, R8309 (1992)] and the Buldyrev et ai. model [Phys. 



Rev. A 45, R8313 (1992)] by Braunstein et aJ. [J. Phys. A 32, 1801 (1999); 
Phys. Rev. E 59, 4243 (1999)]. Even these models have being classified in 
the same universality class that the QKPZ the contributions to the growing 
mechanisms are quite different. The lateral contribution in the DPD models, 
leads to an increasing of the roughness near the criticality while in the QKPZ 



equation this contribution always flattens the roughness. These results sug- 
gest that the QKPZ equation does not describe properly the DPD models 
even when the exponents derived from this equation are similar to the one 

obtained from simulations of these models. 
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I. INTRODUCTION 

The description of the noise-driven growth that leads to self-affine interface far from 
equihbrium is a challenging problem. The interface has been characterized through scaling 
of the interfacial width w with time t and lateral size L. The result is the determination of 
two exponents [3 and a called djTiamical and roughness exponents, respectively. It is well 
known that interfacial width w ^ L"' for t ^ t* and w ^ t^ for t <^ t*, where t* ~ L°/^ 
is the saturation time. These properties occur in many experimental situations and many 
models of surface growth. The values of the exponents leads to their classification in different 
universality classes. Several models, belonging to the same directed percolation depinning 
(DPD) universality class, have been introduced to explain experiments on fluid imbibition in 
porous media, roughening in slow combustion of paper, growth of bacterial colonies, etc. It 
is currently accepted that the quenched disorder plays an essential role in those experiments. 
The DPD models take into account the most important features of the experiments |jl|J^ . On 
the other hand the phenomenological stochastic differential equation proposed by Kardar, 
Parisi, and Zhang with quenched noise (QKPZ) is used to describe the interface growth 
in disordered media and drives to the same universality class than the DPD models, in the 
sense that they have the same exponents. Moreover, processes with the same exponent 
may not belong to the same universality class. For example, 1 + 1-dimensional lattice gas 
simulations of roughening of immiscible fiuid-fiuid interface @] lead to the same exponents 
as the 1 + 1-dimensional KPZ P (/? = 1/3 and a = 1/2) for surface growth, but this model 
is completely linear, so there is no obvious mathematical relationship between these two 
processes. 

The two first DPD models were simultaneously introduced by Buldyrev et al. 0] and 
Tang and Leschhorn |Q to explain the fluid imbibition in paper sheet. Several authors have 
been focused their attention on scaling properties and relationships between the dynamical 
and the static exponent for these models. These two models have been recently reviewed 
by Braunstein et al. p,|^ from a different point of view than the traditional one. The 



principal contribution was the restatement of the Microscopic Equation for each modeL 
These equations allows the separation into two contributions: the local and the lateral (or 
contact) one. They found that the lateral contribution to the temporal derivative of the 
interface width (DSIW) may be either negative or positive and that the behavior of this 
contribution depends on the pressure p, where p is the microscopic driving force. The 
negative contribution tends to smooth out the surface, this case dominate for p ^ Pc (where 
Pc is the critical pressure). The positive contribution enhances the roughness. At the 
critical pressure the local contribution to the DSIW is practically constant, but the lateral 
contribution is very strong. This last contribution, has important duties on the power law 
behavior in the DPD models. 

In this paper we focus the attention in the QKPZ equation in order to compare the 
different contributions with the ones obtained for the DPD models. In this context we show 
that the results obtained from this equation are quite different from the ones obtained in the 
DPD models. We separate the QKPZ equation into two contributions: the local contribution 
and the lateral one. In this context we study the mean height speed (MHS) and the DSIW 
as a function of time. The paper is organized as follows. In Section y we separate the QKPZ 
equation into two contributions for the mean height and the roughness. We study the MHS, 
analyzing the local and the lateral contributions. Also, the two contributions to the DSIW 



are analyzed. In Section III we compare the DPD models with the QKPZ model. Finally, 



we present the main conclusions in Section IV. 



II. MACROSCOPIC CONTRIBUTIONS FROM THE QKPZ EQUATION 

A. Equations 
The QKPZ equation for the surface height h = h{x,t), in 1 + 1-dimension, is given by 

dth = J^ + vdlh + \{d^hf + 7]{x,h) (1) 



where z/ and A are constants and the noise depends on position x and height h with the 
properties that < ri{x, h) >= and < ri{x, h) ri{x', h') >= 2D 6{x — x') 6{h — h') where 2D 
is the noise intensity. 

We can distinguish two contributions to this equation, the local growth S = J-" + ri{x, h) 
and the lateral one C = u d^h + X {dr^hy. So we can write the evolution equation for the 
mean height as 

dth = S + C. (2) 

Taking the derivative of the square interface width, w"^ = {{h— (/i))^), its evolution equation 
is given by 

dtw^ = 2{(h - {h))dth) = 2{(h - {h))S) + 2{(h - {h))C) (3) 

where (...) means average over the lattice. The first term can be identified as the local 
growth contribution, and the second term as the lateral growth contribution. The separation 
into these two analytical terms allows us to compare the mechanisms of growth in the QKPZ 
equation with ones of the DPD models. In the present paper we focus only on the dynamical 
behavior, i.e. t ^ t* ~ L (in these models a ~ /?) for the mean height and roughness. 

We have performed the direct numerical integration of Eq. (|l]) in one dimension in the 
discretized version P,|| 



h{x, t + At) = h{x, t) + At { h{x - 1, t) + h{x + 1, t) - 2h{x, t) 
+ - {h{x + l,t)-h{x-l,t)f 

o 
+J^ + r^{x,[h{x,t)])} , 

where [. . .] denotes the integer part and r] is uniformly distributed in [—a/2, a/2], where 
a = 10^/^ is selected. We use L = 8192 and At = 0.01. The initial condition is h{x, 0) = 
and periodic boundary conditions are used. The averages was taken over 100 samples. We 
study this equation for different values of A with z/ = 1. 



B. Mean height 

In Fig. we show the MHS as a function of time in the three regimes (moving, critical 
and pinning phases) for A = 1. The initial condition for the MHS is T in all regimes. At the 
criticality we found for the mean height a power law behavior with approximately the same 
dynamical exponent that the roughness one, /? = 0.67 ± 0.01 for Tc = 0.464, where J-'c is 
the critical driving force. In the moving and pinning phases we can see that this power law 
does not hold. Bellow the criticality, in the pinning phase, the MHS goes to zero. In the 
moving phase (JF > JF^), the MHS goes to certain constant value. 

In Fig. ^ we show the contributions to the MHS: the local one (S) and the lateral one 
(C). The local contribution, which is equal to JF at t = 0, is stronger in the early time 
regime. This is because in this regime the difference of heights between nearest neighbors 
is very low and the contribution of the lateral term is negligible. We see that the local 
contribution takes negative values from t ^ 7, so in this regime the local contribution brakes 
the growing of (h). The behavior of both contributions, lateral and local one, is not very 
different in every phases although their sum only drive to a power law at the criticality. In 
every phases both contributions are equal at t ~ 2.5. This means that at this time both 
mechanisms are equal independently of JF. Increasing A the lateral contribution is enhanced 
at shorter time, for A = 2 the crossover is at t ~ 1. In the asymptotic dynamic regime both 
contributions are important and neither dominates over the other one. 

C. Roughness 

Fig. ^ shows the temporal DSIW as a function of time for various values of JF and A = 1. 
Here we found that the DSIW increases continuously from zero to a maximum value. The 
power law holds only at the criticality. The DSIW goes asymptotically to zero at the pinning 
and moving phases. The dynamical exponent obtained for the roughness was (3 = 0.66±0.02. 

In Fig. HI, we show the two contributions to the DSIW for different values of JF. The 



local contribution 2{{h — {h))S) to the DSIW is always positive. As JF decreases, this 
contribution also decreases slowly, but always roughen the interface. On the other hand, the 
lateral contribution 2{{h — {h))C) takes negative values in every phases, smoothing out the 
surface. 

In Fig. ^, we plot the DSIW for three values of A at the criticality. The slope /? is 
independent of A although we found some differences in these plots. The scaling dynamical 
regime is reached before at greater values of A. This is due to the fact that the lateral 
contribution, which is the main responsible of the generation of correlations, becomes more 
important earlier for larger values of A. 

III. COMPARISONS WITH THE DPD MODELS 

In this section we present the similarities and differences between the DPD models and 
the QKPZ equation. In previous work Braunstein et al. p,|^ wrote the microscopic equation 
for the TL and Buldyrev models. They identified two separate contributions for the MHS 
and the DSIW: the local and the lateral (contact) one. In the present paper, we have also 
separate the QKPZ equation into these two contributions. 

We found that, in the asymptotic dynamical regime, the behavior of the MHS are similar 
in the QKPZ equation and in the DPD models in every phases @,0]. At short time we found 
a maximum value which depend on A. This maximum value is not found in the DPD 
models. Differences at short time are expected because the QKPZ only describes the scaling 
limit or hydrodynamic limit. The lateral and local contributions to the MHS in the DPD 
models are qualitative and quantitative different from the QKPZ ones. At long time the 
local contribution brakes the growing of (h) in the QKPZ equation (in the DPD models this 
contribution is always positive). Also, the contributions of the QKPZ equation are not so 
much different in each phase, while in the DPD models both contributions are very different 
behavior in each phase. 

In the asymptotic regime, the behavior of the DSIW is similar in the QKPZ equation and 



in the DPD models. In the DPD models p is the initial condition in all regimes, this is due 
to the fact that in the early regime the dynamics is like random deposition with probability 
p [|1^,|11|]. In the QKPZ equation the DSIW increases continuously from zero to a maximum 



value, the macroscopic equation presented by Braunstein et al. p0| , p!2| for the DPD models 
holds in the scaling limit or hydrodynamic limit, but breaks down at short times as was 
expected. The lateral and the local contributions to the DSIW is quite different in each 
model. In the QKPZ the lateral contribution is always negative. This is the main difference 
with the DPD models where for p ^ pc, the lateral contribution is always positive roughening 
the interface PJ^]- ^"^ the DPD models the lateral contribution is enhanced by local growth, 
the lateral growth may also increase the probability of local growth. This crossing interaction 
mechanism makes the lateral growth dominant near the criticality and this is due to the 



fact that the quenched noise is coupled to the dynamic of the interface |T^. In the QKPZ 
model this cross mechanism between contributions is not taken into account because the 
noise is additive. Moreover, the model are said to belong to the same universality class that 
the QKPZ equation because the exponents are quantitatively similar. Even thought the 
behaviors of the contributions to the growth are qualitatively and quantitatively different. 
In the experiments the advancement of the interface is determinated by the coupled effect 
of the random distribution of the capillary sizes, the surface tension and the local properties 
of the flow, so it is not surprising that all these effect give rise to a multiplicative noise. 
This multiplicative noise must be taken into account at the time to pose a model with the 
essential features of the experiment of surface growth in disordered media. In the TL and the 
Buldyrev models the growing rules for the evolution of the local height are strongly coupled 
to the quenched noise in a multiplicative way. In both models the microscopic rules that 
allows the growth from an unblocked cell [^,0 depends in some way on the local slope. In 
that sense this coupled effect is not taken into account in the QKPZ equation. The effect of 
a multiplicative noise as being proposed by Csahok et al. [^ by means of a phenomenological 
equation. They found a crossover between two temporal regimes with (3 = 0.65 to /? = 0.26 
but the value of a ^ 0.47 was obtained over a short range spatial scale. Braunstein et al. 

8 



derived the continuous equation for the TL model [|T^]. This equation as the same term that 
the QKPZ equation but its coefficients depends on the competition between the driving 
force and the quenched noise. In that sense the noise is multiphcative. This results joint to 
the results obtained in this work supports that the QKPZ does not describe fully the DPD 
models even if the exponents are quantitatively similar. 

IV. CONCLUSIONS 

We express the evolution equations of the QKPZ model for the mean height and the 
roughness into two contributions: the local and the lateral one. We found that the contri- 
butions to the growing mechanisms are quite different from the DPD models. In the scaling 
regime, the local contribution to the MHS brakes the growing of {h) in the QKPZ model. 
The lateral contribution to the DSIW is negative in all phases for the QKPZ model smooth- 
ing out the interface. In the DPD models the lateral contribution is always positive for 
P ^ Pc roughening the interface. Nevertheless, the DSIW and MHS gives the same scaling 
exponents in these models, moreover it is not clear why. Our results suggest that the QKPZ 
equation does not describe properly the dynamics of the DPD models even if the exponents 
are similar. 
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FIGURES 
FIG. 1. Plots of T'^ dh/dt vs t for A = 1. The parameter J" is 0.51 (Q), 0.464 (□), and 0.43 
(A). 

FIG. 2. Semi- In plots of the different contributions to the MHS as a function of time for 
different values of J- and A = 1. The circles (Q) represent the local contribution, the triangles 
(A) represent the lateral contribution, and the squares (□) represent the total MHS. The (a) plot 
shows the critical phase J- = 0.464. The (b) plot shows the pinning phase J- = 0.43. The (c) plot 
shows the moving phase J- = 0.51. 

FIG. 3. DSIW as a function of time in the critical, pinning, and moving phases for A = 1. The 
parameter !F is 0.464 (solid line), 0.43 (dashed line), and 0.54 (doted line). 

FIG. 4. Semi-ln plots of the different contributions to the DSIW as a function of time for 
different values of J- and A = 1. The circles (Q) represent the local contribution, the triangles 
(A) represent the lateral contribution, and the squares (n) represent the total DSIW. The (a) plot 
shows the critical phase J- = 0.464. The (b) plot shows the pinning phase J- = 0.43. The (c) plot 
shows the moving phase ^ = 0.54. 

FIG. 5. DSIW as a function of time in the critical regime for different values of A: 0.5 (solid 
line), 1 (doted line), and 2 (dashed line) 
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